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Preliminary  Design  and  Analysis  of  Heading  Control  Systems 
for  a  Jet  Controlled  Test  Vehicle  (RTV1  j) 

by 

A.  P.  MacLaren,  B.Sc. ,  A.M.I.E.E. 
and 

J.  H.  W.  Shannon,  B.E.  (Mech. ),  B.E.  (Aero. ). 


SUMMARY 


As  part  of  the  supporting  experimental  programme  for  ballistic  missile 
development,  it  is  intended  to  control  the  heading  of  uriboosted,  vertically- 
launched,  'wingless  RTV1  missiles  by  means  of  rocket  motor  deflection.  A 
separate  system  of  air  jets  will  be  used  for  roll  stabilisation. 

This  note  describes  in  some  detail  the  preliminary  design  and  analysis 
of  the  proposed  pitch  and  yaw  control  system.  Frequency  response  techniques 
are  used  in  the  basic  design  of  the  system.  Prom  this  analysis,  in  which 
aerodynamic  forces  are  neglected,  suitable  values  of  loop  gain  for  a  range 
of  rocket  motor  servo  bandwidths  are  determined.  Two  possible  systems,  with 
and  without  integral  of  error  control,  are  considered  and  their  relative 
merits  are  illustrated  by  comparing  the  deviations  from  a  vertical  trajectory 
caused  by  an  initial  heading  error  and  a  oontrol  system  bias.  These  simpli¬ 
fied  trajeotory  calculations  indicate  that  the  oontrol  system  including 
integral  of  error  i3  to  be  preferred.  In  the  final  seotion,  the  effeots  of 
aerodynamic  forces,  under  particular  conditions  of  flight,  on  the  stability 
of  the  proposed  system  are  examined. 

It  is  concluded  that  a  stable  heading  control  system  for  aerodynamically 
unstable  RTV1  test  vehicles  can  be  developed  along  the  lines  proposed  in  this 
note,  without  demanding  excessively  high  component  performances.  Further 
work  will  include  studies  of  the  effects  of  body  flexure  on  the  stability  of 
the  system  and  simulator  tests  in  which  some  of  the  actual  components  will 
be  incorporated  in  the  oontrol  loop. 
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1  Introduction 

In  flight  trials  to  study  the  use  of  rocket  motor  deflection  as  a  means 
of  control"' ,  it  is  proposed  that  modified  RTV1  test  vehicles  shall  be 
launched  vertically  from  a  zero  length  launcher.  No  boost  will  be  used  and 
the  missile  longitudinal  acceleration  at  take-off  will  be  less  than  1g. 

For  a  typioal  case  a  thrust  duration  of  22f  seoonds  may  be  expeoted  to  give 
a  speed  of  the  order  of  J00  ft/seo  and  an  altitude  of  about  8,000  ft  at 
all-burnt.  In  these  experiments,  there  will  be  an  initial  period  of  flight 
in  which  the  aerodynamic  forces  will  be  negligible.  In  this  note,  therefore, 
the  analysis  on  which  the  control  system  design  is  based  has  been  simplified 
by  neglecting  aerodynamic  forces  in  the  first  instance.  The  effeot  of 
aerodynamic  forces  on  the  performance  of  the  proposed  control  system,  for  one 
particular  flight  condition,  is  examined  subsequently.  For  these  trials  the 
normal  RTV1  wings  and  fins  will  be  removed  and  consequently  the  modified 
vehicle  will  be  aerodynamically  unstable. 

The  missile  will  be  roll  stabilised  to  a  gyro  datum  by  means  of  a 
separate  system  of  air  jets.  The  design  of  the  latter  is  outside  the  scope 
of  this  note.  For  lateral  control  the  missile  heading  angle  will  be  con¬ 
trolled  in  two  orthogonal  planes,  referred  to  as  the  pitch  ard  yaw  planes. 

As  the  configuration  is  symmetrical,  control  in  one  plane  only  is  considered, 
the  system  in  the  other  plane  being  assumed  to  be  identical.  Cross-coupling 
effects  are  not  included  in  the  analysis  as  no  reliable  data  on  these  are 
available. 

Two  possible  types  of  lateral  control  system  are  considered.  In  the 
first, rocket  motor  deflection  proportional  to  heading  angle  error,  sensed  by 
a  free  gyroscope,  is  used  with  rate  feedback  from  a  rate  gyroscope  for 
damping.  This  is  the  simplest  system  which  will  give  stable  flight.  The 
second  is  similar  to  the  first  with  the  addition  of  a  signal  proportional  to 
the  integral  of  heading  angle  error  to  eliminate  steady  state  errors.  Blook 
diagrams  of  the  two  types  of  system  are  shown  in  Figs.  1  and  2.  frequency 
response  techniques,  which  are  generally  employed  in  the  analysis  of  linear 
servomechanisms,  are  used  in  the  basic  design  of  the  system.  The  transfer 
functions  of  the  rate  gyroscope  and  the  motor  position  control  servo  are 
assumed,  for  this  purpose,  to  be  simple  quadratic  functions.  In  the  design 
of  the  system,  three  possible  motor  servo  bandwidths,  corresponding  to 
natural  frequencies  of  3,  5  and  10  cycles  per  second,  are  examined  to  show 
the  effect  of  this  parameter  on  the  permissible  values  of  loop  gain  for  a 
well  damped  transient  response.  These  bandwidths  are  thought  reasonable  in 
the  light  of  existing  hydraulic  servo  development.  Scan©  relevant  design  work 
is  desoribed  in  Ref.  2. 

In  section  3  the  lateral  motions  arising  from  an  error  in  missile  atti¬ 
tude  at  the  instant  of  launch  and  from  a  constant  rocket  motor  deflection 
bias  are  calculated  for  the  two  types  of  system.  These  calculations  based 
on  certain  simplifying  assumptions,,  indicate  that  the  system  including 
integral  of  error  is  more  effective  in  reduoing  deviations  from  the  desired 
course  and  is  therefore  to  be  preferred.  The  stability  of  thi3  system  with 
missile  aerodynamios  included  is  analysed  in  section  4- 

2  Control  System  Analysis  Neglecting  Aerodynamic  Foroes 
2.1  Response  of  the  missile  to  rooket  motor  defleotion 

Referring  to  Fig. 3,  it  oan  be  seen  that  the  equation  of  missile  nngnin-r 
motion  in  the  yaw  plane  when  aerodynamic  foroes  are  negligible,  may  be 
written 


0?  a  -  Hlb 


(1) 
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where  the  dot  notation  indioates  differentiation  with  respect  to  time 

and  f  <=  the  missile  heading  angle  with  respect  to  a  fixed  direction 
(assumed  vertical), 

5  »  the  rooket  motor  defleotion  measured  from  the  missile  axis 
(assumed  small  so.  that  sin  6  4  6), 

T  a  the  rocket  motor  thrust, 

C  s  the  missile  lateral  moment  of  inertia  about  an  axis  through  the 
centre  of  gravity, 

i  a  the  distance  between  the  missile  centre  of  gravity"  and  the 
rocket  motor  pivot  point. 

Re-writing  equation  (l )  in  the  form  of  a  transfer  function  we  have 


I  » 


(2) 


where 

and  G^ 
and  Z 


b  —  s  the  invariant  factor  of  the  transfer  function,  . 

=  —  b  the  frequency  variant  part  of  the  transfer  function, 

denote  the  Laplace  transforms  of  ljr  and  6  respectively. 


2.2  Transfer  functions  of  the  servo  and  rate  gyroscope 

It  is  assumed  that  the  closed-loop  transfer  function  of  the  rooket  motor 
position  control  servomechanism  oan  be  represented  by  a  quadratic  function,  of 
the  form 


v  2 

6  K3m3 

•ST  =  -  K,G,  a  -  -5 - — - k 

e3  p  ^  p.  +  2£^p  * 


D) 


where  =  the  input  signal  to  the  servo, 

Kj  a  the  invariant  factor  of  the  transfer  function, 

G-j  a  the  frequenoy  variant  part  of  the  -transfer  function. 


2%  x  the  undamped  natural  frequency  f y 
the  damping  factor. 

To  show  the  effeot  of  different  servo  bandwidths  on  the  permissible 
values  of  heading  control  stiffness,  three  values  of  the  undamped  natural 
frequency  f,,  namely  3,  5  and  10  cycles  per  seoond  have  been  used  in  the 
analysis.  The  damping  factor  is  assumed  equal  to  0.6  throughout. 


For  the  rate  gyroscope  a  quadratic  transfer  is  also  assumed,  namely 


K2G2 


V>2 


2,  2. 

p  +  2£2a2p  + 


(4) 
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V 


vj 


=  the  output  signal  from  the  rate  gyrosoope, 

and  Gn  =  the  invariant  and  frequency  variant  parts  of  the 
transfer  function  respectively, 

»  2%  x  the  undamped  natural  frequency  fg, 

=  the  damping  factor. 

The  values  used  for  fg  and  are  respectively  29  cyoles  per  second 
and  0,5.  These  are  related  to  an  existing  rate  gyrosoope  which  has  a  range 
of  ±60°  per  second. 

Fig, 4  shows  the  frequency  response  curves  corresponding  to  the  transfer 
functions  assumed  for  the  servo  and  rats  gyro.  The  values  of  the  gain 
parameters  Kg  and  Kj  are  assumed  to  be  at  our  disposal. 

2. 3  Frequency  response  study  of  oontrol  system  without  integral  of  error 

It  oaii  be  seen  from  the  block  diagram  of  Fig.1  that  the  inner  loop 
feedback  transfer  function  is  given  by 


where  e2 
*2 

“2 

% 


h  J 

¥  S 


JL 

53 


*  *2V-  K3G3 


V4 


(5) 


The  frequency  variant  part  (GgG^G^  the  transfer  function  of  expression 
(5)  is  plotted  in  Fig.  5.  Three  3ets  of  curves  are  shown  corresponding  to 
the  three  assumed  servo  frequencies  of  3>  5  and  10  cycles  per  second.  It 
can  be  seen  by  plotting  gain/phase  curves  on' a  Nichol's  chart  that  the 
closed  loop  response  of  the  rate  feedback  loop  will  be  adequately  damped  if 
a  phase  margin  of  60°  is  maintained  in  each  case.  For.  this  condition  the 
following  loop  gains  and  stability  margins  are  obtained. 

Table  I 


Rate  Feedback  Loop  Gains 


Servo  bandwidth 
f-j  cycles  per  second 

3 

5 

10 

Loop  gain  in  db  (20  log^^L^) 

17.0 

21.0 

26.0 

Loop  gain  factor  L^  =  K^K^K^ 

7.08 

11.2 

.  20.0 

Loop  gain  margin 

9.1  db  at 
2.82  o/s 

8. 6  db  at 
4-52  o/s 

8.3  db  at 
8.25  o/s 

Loop  phase  margin 

60°  at 
1.14  o/s 

60°  at 

1.81  c/s 

60°  at 

3.26  o/s 

With  the  rate  feedback  loop  closed  we  have 


and  the  closed  loop  transfer  function  may  be  expressed  as 


SECRET 


7 


(6) 


SECRET 


Technioal  Note  No.  GW  39 7 


i  1  r  wwt  ~ 
s1  =  ^2  L> +  w3  w J 


where  a  the  output  signal  from  the  heading  angle  gyroscope.  The  feedback 
transfer  function  for  the  position  feedback  loop  is  given  by 


I 

*S 


K. 


pK„G, 


-  WSWk 


'2  2  L' 


1  + 


(8) 


where  i|r  denotes  the  heading  angle  error  from  the  gyro  datum  and  Hj  is  the 
gain  factor  for  heading  angle  error  signal. 

Gain  and  phase  plots  representing  the  frequency  variant  part  of  the  above 
transfer  function,  incorporating  the  appropriate  rate  loop  gains  given  in 
Table  I  for  the  three  cases,  are  shown  in  Pig.  6.  Inspection  of  these  on  a 
Nicholas  chart  shows  that  a  well  damped  closed  loop  response  will  be  obtained 
with  a  phase  margin  of  60°.  The  following  loop  gains  and  stability  margins 
are  obtained  for  this  condition. 


Table  II 


Position  Feedback  Loop  Gains 


# 


u 


* 


Servo  bandwidth 
f ^  cycles  per  second 

3 

5 

10 

Loop  gain  in  db  (20  log^^L^) 

11.0 

14.8 

19.6 

K1 

Loop  gain  factor  1^  =  ^- 

3.55 

5.50 

9.55 

Loop  gain  margin 

8.7  db  at 
1.60  c/s 

9.1  db  at 
2.48  c/s 

8.9  <3b  at 
4.30  c/s 

Loop  phase  margin. 

60°  at 
0.58  c/s 

60°  at 
0.88  c/s 

60°  at 
1,50  o/s 

Pig.  7  shows  the  variation  with  frequency  of  the  gain  of  the  closed  loop 

f 

transfer  function  -ijr-  for  the  three  cases.  ^  denotes  the  gyro  datum  angle 

measured  with  respect  to  the  same  fixed  direction  as  ^r, 

2,4  Frequency  response  study  of  control  system  with  integral  of  error 

As  the  rate  feedback  loop  transfer  function  is  identical  with  that  already 
studied  in  seotion  2.3,  the  loop  gains  set  out  in  Table  I  also  apply  to  the 
control  system  illustrated  by  Pig.  2.  In  this  oase,  however,  the  position 
feedback  loop  includes  an  additional  term  proportional  to  integral  of  heading 
angle  error  and  hence  the  feedback  transfer  function  of  this  loop  beoomes 
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*2 


/  K1' 


2„ 

P  Go 


:)r  WWi 

U  +  WjWi 


(9) 


where  denotes  the  gain  factor  of  the  integrator. 


The  frequency  response  of  the  term 


P\  L1  +  wywu 


is  shown 


in  Pig.  8.  The  parameters  of  the  remaining  part  of  the  transfer  function  of 

K. 

expression  (9)  (l  +  P  v~  )  axe  our  disposal  to  be  chosen  in  such  a  way 


K1+*4) 


that  the  response  of  the  outer  loop  is  satisfactory.  Sinoe  there  are  two 
parameters  at  our  disposal,  namely  K-j  and  Kc,  the  transfer  function  of 
expression  (9)  may  be  made  to  satisfy  two  conditions.  To  facilitate  com¬ 
parison  between  the  types  of  system  shown  in  Pigs.  1  and  2,  the  criterion  of 
60°  phase  margin  in  the  outer  loop,  used  in  section  2,3,  is  retained  as  the 
first-  of  these  oorxLitions.  The  second  condition  selected  is  that  Kc,  the 
integrator  gain  faotor,  shall  be  as  large  as  possible,  compatible  with  the 
first  condition  being  fulfilled,  to  reduce  errors  due  to  low  frequency  dis¬ 
turbances  to  a  minimum. 

The  method  used  to  determine  the  values  of  K5  and  JCj  to  satisfy 
these  conditions  is  as  follows.  Consider  first  the  frequency  response  of 

/ 

the  term  (1  +  p  •=-  )  illustrated  in  Pig.  9*  This  has  a  break  frequency,  as 
yet  undetermined  at  <y  a  <=*■  .  Above  this  frequency  the  gain  asymptote 

Q  ^ 

rises  at  6  db  per  octave  and  the  phase  lead  increases  from  45°  to  9O0.  Let 
and  denote  respectively  the  gain  (in  db)  and  phase  of  the  term 
K.j  , 

at  the  outer  loop  cut-off  frequency  uQ  i-e.  the  frequency  at 

which  the  feedback  gain  is  unily  (zero  db).  Then  to  fulfil  the  60°  phase 
margin  condition  the  following  equations  must  be  satisfied:- 


20  log10  +  D1  +  D2  *  0  (db) 
+  j a  ~  120  (degrees) 


(10) 

(11) 


where  Dg  and  denote  respectively  the  gain  in  db  and  the  phase  of 
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at  the  out-off  frequency  « 


The  out-off  frequency  Uq  corresponding  to  any  particular  value  of  jij 
is  easily  determined  by  means  of  equation  (11)  and  the  phase  curve  of  Pig.  8. 
The  gain  D2  corresponding  to  any  value  of  oan  be  read  from  the  gain 

curve  of  Pig.  8.  For  example,  for  ^  a  80°,  a  -200°  from  equation  (11) 

and  from  Pig.  8  it  oan  be  seen  that  Uq  a  3*72  radians  per  second  and  D2  a 
-22.7  db  for  a  servo  natural  frequency  of  5  pyoles  per  seoond.  The  nine 
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values  of  wc  aid  the  corresponding  values  of  Dg  for  the  three  servo  band- 
widths  and  for  <f>\  —  75° >  80°  and  83°  are  given  in  Table  III,  Prom  Pig.  9  it 
can  be  seen  that  the  values  of  Dj  corresponding  to  ^  »  75°,  80°  and.  85° 
are  1.1.8  db  at  a  frequency,  of  3*7 3  Wq>  15«3  db  at  5*8  w0  and  21.2  db  at 
11.6  wQ.  Since  these  are  the  cut-off  frequencies  we  have:- 


w  w  w„ 

.  c  c  _  c 


=  wQ  a  -jTyj  ,  Jjr  or  for  *  75  ,  80  or  85  respective^. 

1  02) 

*5 

The  nine  values  of  the  ratio  given  in  Table  III  are  derived  from 

*1 

equation  (12)  by  using  the  appropriate  value  of  to  and  the  corresponding 
Kc  ° 

values  of  rf-  are  obtained  from  equation  (lO).  Hence  the  values  of  K/t  and 

k2 

Kr  corresponding  to  the  three  values  of  j>^  are  determined  since  Kg  is 
fixed  by  the  rate  feedback  loop  gain  given  in  Table  I. 

Table  III 

Values  of  and  corresponding  to  the  three  values  of  ^ 


Phase  ^  degrees 


Gain  Oj  db 


Wo 

Ratio  — 
W„ 


Phase  degrees 


3.75 

-195 


5.80 

-200 


Servo  bandwidth 
r3  c/s 

3 

5 

10 

3 

5 

10 

3 

5 

10 

Cut-off  frequency 
ej,  rad/sec 

1.80 

2.85 

4.85 

2.33 

3*72 

6.35 

2.90 

4.65 

8.00 

to.  (-  -—l  rad/sac 

w 

Gain  0,  db 

0.480 

0.760 

1.29 

0.402 

0.642 

1.10 

0.250 

0.402 

0.690 

-10.2 

-18.2 

-27 .3 

-14.7 

-22.7 

-32.2 

-18.5 

-26.5 

-36.2 

k5 

20  lo&jQ  ^ 

-1.6 

*6.4 

♦15.5 

-0.6 

♦7.4 

♦16.9 

-2.7 

♦5.3 

♦15.0 

Kc 

Ratio  -2 

k2 

0.832 

2.09 

5.96 

0.933 

2.J4 

7.00 

0.733 

1.84 

5.62 

Ratio  51 
*2 

1.73 

2*75 

4.62 

2.32 

3.64 

6.37 

2.94 

4.58 

8.14 

11.6 

-C05 


*5  ^ 

All  the  values  of  rf-  and  tr-  given  in  Table  III  satisfy  the  first 
*2  *2 

(60°  phase  margin)  condition.  Comparison  of  the  values  of  =-  ,  given  in 

*2 

Table  II,  for  the  system  without  integral  of  error  with  the  corresponding 
values  in  Table  III  shows  that  the  latter  are  lower.  For  predetermined  values 
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of  K £,  IC)  increases  with  fa  tending  to  the  values  given  in  Table  II  as 
fa  tends  to  J0°  and  Kc  tends  to  zero.  It  can  be  seen  from  Table  III, 
however,  that  for  all  three  servo  bandwidths  is  larger  for  fa  =  80° 
than  for  fa  =  75°  or  85°.  Thus  has  a  maximum  value  in  the  region  of 
fa  =  80°.  The  values  of  Kj  and  Kc  given  in  the  three  central  columns 
therefore  satisfy  approximately  both  the  chosen  conditions  and  will  be  used 
in  all  subsequent  calculations  on  the  integral, of  error  system. 

The  process  of  determining  the  values  of  Kc  and  Kj  which  satisfy 

liclay 


the  conditions  selected  can  be  carried  out'  quickly  and  conveniently  by 
using  templates  cut  to  the  shape  of  the  gain  and  phase  curves  of  Fig.  $. 

The  frequency  response  curves  for  the  outer  loop  feedback  transfer 
function  for  the  three  servos  and  for  fa  =  80°  are  shovm  in  Fig.  10.  It 
can  be  seen  from  this  that  the  gain  margin  is  12.4  db  at  1.56,  2.42  and 
4.15  cycles  per  second,  the  phase  margin  being  60°  at  0.37>  0*59  and  1.01 
cycles  per  second,  for  f  j  =  J>,  5  and  10  cycles  per  seoond  respectively. 

Fig.  11  shows  a  comparison  of  the  frequency  responses  for  fa  =  80°  and 
fa  =  85°  for  a  10  cycles  per  second.  It  is  evident  that  fa  s  80° 
gives  a  higher  gain  at  low  frequencies,  due  to  the  higher  integrator  gain 
factor,  at  the  expense  of  a  rather  lower  gain  at  high  frequencies.  The 
variation  with  frequency  of  the  gain  of  the  closed  loop  transfer  funotion 

f 

for  the  three  servo  bandwidths  is  shown  in  Fig.  12.  Polar  plots  for 

V 

f^  n  10  cyoles  per  seoond  for  the  two  types  of  system  shown  in  Figs.  1  and  2 
are  compared  in  Fig.13. 

2*5  Numerical  values  of  gain  factors  for  the  proposed  control  systems 

The  values  given  in  this  section  ere  based  on  the  frequency  response 
studies  of  sections  2.3  and  2.4  and  on  the  following  data  for  the  RTV1  missile. 

Table  IV 

Assumed  Missile  Parameters 


Fuel  Mass 

Condition  m  slugs 

Moment  of 
Inertia 

C  slug  ft2 

Centre  of 
gravity 
position  • 
aft  of 
nose 
ft 

Rocket 
motor 
pivot  aft 
of  centre 
of  gravity 
l  ft 

Assumed 
Motor 
Thrust 
T  lb 

Tanks  full  1 7. 6 

352 

8.94 

7.79 

1000 

Half  full  15.6 

346 

8.85 

7.88 

1000 

Tanks  empiy  1 3*5 

328 

8.55 

8.19 

1000 

h  °-o 

radians 
per  s  sem¬ 
per  rad 

radians 
per  sec^ 
per  deg 

22.1 

22.8 

25.0 

0.386 

0.397 

0.436 

In  calculating  the  values  of  L,  Kg  and  K^  given  in  the  following 
tables,  Kj  is  assumed  equal  to  uniTy  (one  degree^of  motor  deflection  per 
volt).  The  actual  value  assumed  is  unimportant  since,  for  given  values  of 
and  the  feedback  loop  gains,  the  products  Kj  Kj,  Kg  K3  and  Kt;  K*  will 
be  constant.  By  choosing  Kj  as  unity,  the  value,  for  ins  tan  oe,  of  Kj  in 
volts  per  degree  is  numerically  the  same  as  that  of  the  product  Kj  Kj  in 
degrees  of  motor  deflection  per  degree  of  heading  angle  error. 

The  parameters  for  the  control  system  of  Fig.1  are  based  on  the  loop 
gains  given  in  Tables  I'  and  II. 
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Table  V 

Parameters  for  the  Control  System  of  Pig. 1 


Fuel  ' 
Condition 

*3 

deg/volt 

*4 

sec“2 

fj  =  3  cycles 
per  seo 

f 2  =  5  cycles 
per  seo 

f-,  a  10  cyoles 
"'per  seo 

*2 

V/deg/sec 

V/deg 

V/deg/sec 

V/deg 

*2 

V/deg/seo 

7 

V/deg 

Pull 

1.0 

22.1 

0.318 

1.13 

0.507 

2.79 

0.905 

8.64 

Hall  full 

1.0 

22.8 

0.311 

1.10 

0.493 

2.72 

0.880 

8.40 

Empty 

1.0 

25.0 

0.281k 

1.01 

0.449 

2.47 

0.800 

7.64 

The  parameters  for  the  control  system  of  Pig.  2  are  based  on  the  gain  ratios 
R,-  K. 

and  given  in  the  central  columns  of  Table  III.  The  values  of  K?,  Kit 
*2  *2  17 

and  are  the  same  as  those  given  in  Table  V  above. 

Table  VI 

Parameters  for  the  Control  System  of  Pig.  2 


Fuel 

Condition 

f  -1  -  3  qyoles 
per  seo 

f  j  =  5  qyoles 
per  sec 

fj  as  10  cycles 
per  sec 

■  h 

.V/sec/deg 

*» 

V/deg 

*5 

V/sec/deg 

*1 

V/deg 

*5 

V/sec/deg 

.*1 

V/deg 

Full 

0.297 

0.738 

1.19 

1.84 

6.33 

5.76 

Half  full 

0.290 

0.722 

1,15 

1.79 

6.16 

5.60 

Empty 

0.265 

0.658 

1.05 

I.64 

5.60 

5.10 

It  can  be  3een  from  the  above  tables  that,  if  the  feedback  gain  factors  are 
kept  constant  throughout  flight,  the  loop  gains  will  increase  as  the  flight 
proceeds  by  approximately  1j$,  i. e.  about  1  db,  due  to  reduction  in  the  moment 
of  inertia  etc.  This  is  unlikely  to  have  serious  effects  in  view  of  the  large 
gain  and  phase  margins  available.  In  the  following  sections  of  this  note, 
therefore,  it  is  assumed  that  the  feedback  gain  factors  Kj,  Ify,  K*  and  will 
be  kept  oonstant  at  the  values  given  in  Tables  V  and  VI  for  tanks  full 
conditions. 

3  Effect  of  Control  System  Parameters  on  Test  Vehicle  Behaviour 

The  optimum  control  system  may  be  defined  as  that  which  allows  the  missile 
to  deviate,  under  the  influence  of  external  disturbing  forces  and  imperfections 
in  the  missile,  by  the  least  amount  from  the  desired  trajectory.  Calculation  of 
trajectories  resulting  from  a  variety,  of  input  demands,  initial  conditions, 
disturbing  forces,  winds,  etc,  is  a  formidable  task  which  may  be  more  con¬ 
veniently  performed  by  means  of  electronic  oomputors.  In  this  work,  therefore, 
no  attempt  has  been  made  to  determine  the  optimum  control  system  as  defined 
above,  but,  in  order  to  illustrate  the  effeot  of  the  oontrol  system  parameters 
on  the  trajectory,  the  missile  motion  resulting  from  two  types  of  imperfection 
has  been  calculated. 

These  illustrative  calculations  are  based  on  the  system  parameters  given  in 
Tables  IV,  V  and  VI  of  section  2.5  and.  on  the  following  simplifying  assumptions 
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(i) 

that  the  demanded  heading  angle  i|k  remains  vertioal  throughout 
flight ; 

* 

(ii) 

that  the  angles  *  and  8  are  small.  The  rocket  motor  deflea 
tion  8  is  assumed  to  be  limited  to  i5°» 

■* 

(iii) 

that  the  missile  mass,  centre  of  gravity,  moment  of  inertia  and 
feedback  gain  factors  remain  constant  at  the  value  corresponding 
to  tanks  full  conditions; 

(iv) 

that  aerodynamio  forces  can  be  neglected; 

•  (v) 

that  the  rate  gyroscope  and  rocket  motor  servo  lags  can  be 
neglected. 

The  deviation  of  the  missile  from  the  Heal  vertical  trajectory  has 
been  calculated  for  the  following  two  oases:- 

(a)  the  missile  ha.s  an  initial  heading  angle  error  of  3°  with 
respect  to  the  vertioal  datum ; 

(b)  a  constant  rocket  motor  deflection  bias  of  1°  is  present  in  the 
control  system. 

The  details  of  these  calculations  are  set  out  in  Appendix  I  and  the 
w  results  are  summarised  in  seotions  3*2  and  3*3* 

3*1  Basic  equations 

’  The  assumptions  stated  above  enable  the  missile  motion  to  be  expressed 

in  terms  of  low-order  linear  differential  equations  with  constant  coeffi¬ 
cients,  the  solution  of  which  is  readily  obtainable  by  the  Laplaoe  Transforma¬ 
tion  method.  Erom  equation  (1 ),  the  angular  motion  neglecting  aerodynamic 
forces,  is  given  by 


V  ~  6  -  -  6  0?) 

where  is  constant  and  equal  to  22,13  by  assumption  (iii). 

By  assumption  (v),  the  control  equation  for  the  system  without  integral 
of  error  (Pig. 1 )  can  be  written 

6  *  -  K3  (e1  -  e2)  n  -  K1K3+e  +  *  (14) 

Now 

83  ^  05) 

since  by  assumption  (i)  ^  a  0,  both  and  ijj)  being  measured  with 

respect  to  a  vertical  datum.  Hence  (14)  becomes 

6  »  +  lyLjjr  .  (16) 

Similarly  the  control  equation  for  the  system  with  integral  of  error 
reduces  to 
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6 


¥3  / 


♦  at  +  KjKjf  + 


(17) 


The  horizontal  component  of  missile  acceleration  is  given  by 

x  =  |  sin  (i)r+  6)  *■“(♦+ 5)  (18) 

by  assumption  ( ii) ,  and  the  vertioal  component  by 

y  *»'  f  COS  (t+  6)  -  g  *  (m  “  g)  *  09) 

The  coordinates  (x,y)  denote  the  position  in  space  of  the  missile  centre  of 
gravity  with  respect  to  its  position  at  the  instant  of  launch;  g  is  the 
acceleration  due  to  gravity  and  m  the  mass  of  the  rocket. 

3 «-2  Missile  deviation  due  to  a  3°  initial  heading  angle  error 

(a)  Control  system  of  Fig.  1  (without  integral  of  error) 

The  results  of  the  trajectory  calculations  for  this  oase,  shown  graphically 
in  Pig.  11,  are  given  in  Tables  VIII  and  IX  of  the  Appendix.  It  can  be  seen  from 
the  latter  that,  although  the  initial  heading  error  decreases  exponentially  to 
zero,  the  transient  angular  motion  results  in  a  steady  state  crossing  velocity 
i.  The  heading  angle  error  is  reduced  to  of  its  initial  value  in  1.20,  0.777 
and  0.526  sec  for  the  loop  gains  chosen  for  f  1  a  3,  5  and  10  c/s  respectively 
and  the  corresponding  steady  state  crossing  velocities  are  0.834,  0.573  and 
0.2*86  ft/sec.  If  the  motor  deflection  limit  is  not  reached,  the  crossing  velo¬ 
city  is  proportional  to  the  initial  heading  angle  error. 

(b)  Control  system  of  Pig. 2  (with  integral,  of  error) 

The  results  for  this  oase  are  illustrated  in  Pig.  15  and  given  in  numerical 
foim  in  Tables  X  and  XI  of  the  Appendix.  Again  the  heading  error  decreases 
exponentially  to  zero  but  the  response  of  this  system  is  slower.  Por  instance, 
the  heading  angle  error  is  reduced  to  2/  of  its  initial  value  in  5»55>  3*48  and 
1.75  sec  for  fj  =  3,  5  and  10  c/s  respectively.  On  the  other  hand,  unless  the 
motor  deflection  limit  is  reaohed  during  the  initial  transient,  no  steady  state 
crossing  velocity  x  occurs  with  this  system.  A  small  position  error  arises 
instead.  The  crossing  rates  of  0.02*7  and.  0.2*6  ft/seo,  which  occur  with  the  gain 
parameters  corresponding  to  f 3  =  5  and  10  o/s  respectively,  are  due  to  the 
effects  of  limiting  the  motor  deflection. 

3»3  Missile  deviation  due  to  a  1°  motor  defleotion  bias 


(a)  Control  system  of  Pig.i  (without  integral  of  error) 

The  trajectories  calculated  for  this  case  are  illustrated  in  Pig.l6  and 
given  in  functional  form  in  Table  XII  in  the  Appendix.  It  can  be  seen  that 
with  this  system  a  1°  bias  causes  steady  state  heading  errors  of  0,88°,  0 .36° 
and  0.11 6°  for  the  gain  parameters  corresponding  to  tj,  a  3,  5  aixL  10  o/s 
respectively.  There  is,  in  consequence,  a  steady  transverse  acceleration  2 
of  0.87,  0.353  and  0.115  ft/seo^.  The  magnitude  of  this  acceleration  is 
directly  proportional  to  the  bias  angle. 

(b)  Control  system  of  Pig, 2  (with  integral  of  error) 

The  trajectories  for  this  oase  are  shown  in  Pig.  17  and  given  in  Table  XIII 
in  the  Appendix.  It  can  be  seen  from  the  form  of  expression  (1.31)  that  no 


-14- 

SECRET 


SECRET 


Technical  Note  No.  GW  397 


steady  state  heading  angle  error  results  from  a  oonstant  bias  in  this  type 
of  control  system.  Consequently  there  is  no  steady  state  transverse  accel¬ 
eration  x  but,  as  can  be  seen  from  Table  XIII  the  transient  angular  motion 
results  in  a  steady  crossing  velocity  t  of  3«30»  0,83  and  0,137  ft/sec  for 
the  gain  parameters  corresponding  to  f 3  =  3,  5  and  10  c/s  respectively. 

This  crossing  rate  is  directly  proportional  to  the  bias,  provided  that  limi¬ 
ting  does  not  occur. 

3.4  Discussion  of  results 

It  has  been  shown  that  without  an  integral  of  heading  angle  error  term 
in  the  control  equation  a  steady  state  component  of  velocity  normal  to’  the 
desired  trajectory  arises  from  an  error  in  launching  attitude,  whereas  with 
integral  of  error  control  a  position  error  arises  but  no  steady  crossing 
velocity,  provided  that  motor  deflection  limiting  does  not  occur.  Similarly 
a  bias  in  the  control  system  causes  a  steady  state  transverse  acceleration 
without  integral  control  whereas  with  integral  control  a  transverse  velocity 
arises  but  no  steady  lateral  acceleration.  Thus  it  appears  that  the  system 
with  integral  of  error  control,  in  spite  of  its  slower  transient  response 
characteristics,  tends  to  give  smaller  deviations  from  the  desired  trajectory 
due  to  disturbances.  The  advantage  of  the  integral  oontrol  system  becomes 
more  marked  as  the  time  of  controlled  flight  increases  since  the  ratio  of  the 
lateral  displacements  arising  from  a  given  disturbance  increases  with  time. 

The  effect  of  the  higher  loop  gains,  which  are  permissible  with  increased 
servo  bandwidth,  can  be  seen  from  Pigs.  14,  13,  16  and  17»  In  all  cases, 
where  limiting  does  not  seriously  affect  the  result,  the  displacement  from 
the  desired  trajectory  due  to  a  given  disturbance  becomes  less  as  the  band¬ 
width  and  gain  increase.  Pig, 15  shows  an  example  (fj  »  10  o/s)  of  how  this 
general  tendency  can  be  offset  by  the  adverse  effeots  of  motor  deflection 
limiting.  It  appears  desirable,  therefore,  to  avoid  using  such  high,  values 
of  loop  gain  that  the  probable  range  of  disturbances  will  cause  excessive 
limiting. 

4  Control  System  Analysis  Including  Aerodynamic  Forces 

The  following  analysis  is  intended  to  illustrate  the  effeots  of  aero¬ 
dynamic  forces  awl  momsnts  on  the  stability  of  the  integral  of  error  oontrol 
system.  Since  the  vehicle  is  symmetrical  and  roll  stabilised,  the  analysis 
may  be  simplified,  as  in  section  2,  by  considering  motion  in  one  plane  only. 
Aerodynamic  and  jet  oontrol  f oroes  are  expressed  with  reference  to  missile 
axe 3,  gravitational  forces  and  effects  arising  from  varying  flight  path  velo¬ 
city  being  neglected.  Theoretical  estimates,  by  Chapman,  of  the  normal  force 
coefficient  and  centre  of  pressure  position  for  the  wingless  RTV1  configura¬ 
tion  at  subsonic  speeds  show  non-linear  variations  with  incidence  of  the 
normal  force  and  yawing  noment.  For  the  purpose  of  the  following  dynamio 
stability  analysis,  we  have  considered  the  effects  of  small  perturbations 
about  fixed  values  of  inoidenoe  and  have  assumed  linear  relations  over  this 
limited  range.  The  missile  is  aerodynamically  unstable  for  small  angles  of 
incidence. 


4*1  Response  of  the  missile  to  rooket  motor  def leotlon 
The  basic  equations  of  motion  in  the  yaw  plane  axe 


m  (v  +  Ur) 
Cr 


(2Q) 


where  v  awl  r  denote  the  sideslip  velocity  and  rate  of  turn  respectively, 
m  the  mass  and  U  the  flight  path  velocity.  X  and  N  denote  respectively 
the  normal  force  and  moment  acting  in  the  yaw  plane.  Taking  total  differen¬ 
tials  of  equations  (20)  gives 
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m  (dv  +  Udr)  «  fp  <iv  +  -g  dr  +  ^  d6 


°<*  -  f  ^v+f  dx+f  d6. 


aN 


(21) 


vy- 

Writing  «  Y^,  -gg  =  Xg  etc,  negleoting  Y  and  Ny  and  denoting 
perturbations  by  dashes,  i.e.  dr  a  r',  dv  a  v'  eto,  equations  (21)  become 


m  (v*  +  Ur')  =  Yvv'  +  Y&8' 


Of* 


N 


.V'  t  if.O- 

y*  +  N-6'.  J 


(22) 


Transforming  equations  (22)  and  eliminating  v*  gives  the  transfer 
function 


-»  1 

£-  »  ±- 
-I  e 

6  S 


«  "  _2 


nSp  +  Vs  ~  ?yn6 
p2.-  y^P  +  Uny 


(23) 


N 


where  y  a  — 
v  m 


N. 


»  % 


m 


n  c 


and  r  =  ^ . 


The  estimated  values  of  the  coefficients  of  expression  (23)  are  given  in 
Table  VII  for  U  *  300  ft/sec  and  incidences  of  0°  and  10°.  Missile  conditions 
assumed,  corresponding  to  tanks  half  full,  are  m  a- 15.6  slugs,  C  «  346  slugftr 
and  centre  of  gravity  8*83  ft  aft  of  the  no3e  (Table  IV),  since  the  speed  half 
way  through  powered  flight  will  be  in  the  region  of  300  ft/seo. 

Table  VII 


Aerodynamic  Derivatives  for  Wingless  RTV1  at  300  ft/seo 


Incidence  £ 

yr 

nv 

y5 

n6 

0° 

10° 

-0.0185 

-0.113 

-0.0138 

-0.0117 

64.2 

64.2 

-22.8 

-22.8 

.  Substitution  of  these  values  in  (23)  gives  for  (3  a  0° 


*  ■  fefo  f  1.307)  B  -  0.317  (1  ±  17.45P) 

6  P2  +  0.O183P  -  4.14  (0.494P  -  1)(0*489p  +  t) 

and  for  0  a  10° 


£.  cs  3-33) 

6  P2  +  0.113P  -  3.51 


-  0.95  (1  ±  6«84p)  # 

(0*53Cp  -  l)(0*518p  +  1) 


(24) 
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4 


V 


v 


* 


Denoting  the  transfer  function  -^r  ,  by  -K^G^  we  have  a  0.317 

5 

for  zero  incidence  and  0.95  for  10°  incidence.  Equations  (24)  and  (25)  give 
the  factors  of  the  function  G^  for  the  two  cases  and  Pig.  18  shows  the 
corresponding  frequency  response  curves. 

4*2  frequency  response  study  of  control  system  with  integral  of  error 

Considering  first  the  rate  feedback  loop,  it  can  be  seen  from  Pig. 2 
that  the  open  loop  transfer  function  is  given  by 


(26) 


The  dashes  indicate  perturbations  about  steady  values.  Pig. 19  shows  the 
frequency  response  curves  for  % GoKgGg  for  a  servo  bandwidth  of  5  o/s  and 
the  corresponding  value  of  K2K3  (=  O.507)  given  for  tanks  full  conditions 
in  Table  V  of  seotion  2,  Only  the  5  c/s  case  is  examined  in  detail  here  to 
avoid  complicating  the  diagrams.  Pig. 20  gives  the  frequency  response  curves 
corresponding  to  equation  (26)  for  mean  incidence  values  of  0°  and  10°  and 
Pig. 21  shows  the  corresponding  inner  loop  Nyquist  diagrams. 


-1 

62 


V 


v 


W5W4 


It  can  be  seen  from  Pig.  20  that  the  phase  margin  is  61°  at  1.86  c/s 
for  both  values  of  incidence  and  the  gain  margin  8.6  db  at  4*56  c/s. 
Comparison  with  the  corresponding  figures  for  no  aerodynamics  given  in 
Table  I  shows  that  the  aerodynamio  forces  and  moments  have  a  negligible 
effect  on  the  stability  of  the  rate  feedback  loop.  The  transfer  function 
^G^gG/K^G^  has  a  single  pole  in  the  right  hand  half  of  the  p-plane 
arising  from  the  unstable  aerodynamic  transfer  function.  Since  the  Nyquist 
diagram  does  not  encircle  the  point  (-1  +  jO)  the  divergent  instability  is 
unaltered  by  the  rate  feedback  loop  and  consequently  the  olosed  loop  transfer 
function  has  a  single  pole  in  the  right  hand  half  of  the  p-plane. 


The  feedback  transfer  function  of  the  outer  loop  is  given  by 


31 


*>  I 

r  m  -1 

2 

P 

L1  +  w 

(27) 


Pig.  22  shows  the  frequency  response  curves  of  the  function  in  square 
braokets,  obtained  from  Pig.  20  with  the  aid  of  a  Nichols  chart,  and  Pig.  23 
gives  the  curves  corresponding  to  equation  (27)  for  Kj  »  1.84  and  Kg  =  1.1 9 
(see  Table  Vi).  It  oan  be  seen  from  the  latter  figure  that  the  phase  margin 
i3  57°  at  zero  incidence  and  59°  at  10°  incidence  at  0 .56  cycles  per  second, 
compared  with  60°  at  0.59  c/s  in  the  absence  of  aerodynamic  forces.  The 
gain  margin  is  13*7  <3b  at  2.66  o/s  compared  with  12.4  db  at  2.4 2  c/s  in  the 
absence  of  aerodynamic  foroes. 

The  outer  loop  Nyquist  diagram  is  shown  in  Fig. 24.  It  is  evident  from 
the  complete  diagram  (inset)  that  the  contour  makes  one  anti-clockwise  and 
two  olookwise  (negative)  encirolements  of  the  point  (-1  +  50)  giving  a  net 
total  of  one  negative  enoirclement  (N  a  -1 ).  Since  the  open  loop  transfer 
function  has  a  single  pole  in  the  right  hand  half  of  the  p-plane  (P  a  1)  wa 
have 


Z 


P  +  N 

0 


(28) 
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where  Z  denotes  the  number  of  zeros  of  the  denominator  of  the  closed  loop 
transfer  function.  Hence  the  system  is  stable  with  the  proposed  values  of  loop 
gains. 

It  should  be  noted  that  if  the  outer  loop  gain  were  reduced  by  factors  of 
9.4  (l9«5  db),  at  zero  incidence,  and  12.1  (21.7  db),  at  10°  incidence,  the 
system  would  beoome  unstable  sinoe  there  would  then  be  one  positive  encircle¬ 
ment  of  the  oritioal  point  and  consequently  Z  would  be  equal  to  2.  Thus  the 
system  is  conditionally  stable  in  the  presence  of  aerodynamic  f oroes  -but,  for 
the  two  oases  considered,  the  gain  margins  in  both  directions  are  large. 

4.3  Response  to  a  step  function  change  of  demanded  heading  angle 

It  is  of  interest  to  compare  the  response,  with  and  without  aerodynamic 
forces,  of  the  proposed  system  to  a  step  ohange  of  heading  angle  demand.  The 
calculations  are  simplified  by  neglecting,  in  both  cases,  the  servo  and  rate 
gyroscope  lags. 

In  the  absence  of  aerodynamic  forces  the  response  of  the  integral  of  error 
system  to  a  step  function  of  magnitude  ijjj  is  given  by 


(L,L2p  *  Lfy) 

3  1  1  2,  ~ ' 

P  (p  +  +  L^p 


(29) 


where  L.  a 


'1  -VA>  h-t; 


and  L. 


3 


*2  ’ 


It  is  assumed  that  all  initial  conditions  are  zero  and  that  no ,  limiting 
occurs.  Similarly  in  the  presence  of  aerodynamic  forces 


V 


(L^p  +  L1L6)(n5p  +  iyr6  -  y^)  jfp 


P  [\p(p  -  y^p  +  Uny)  +  (n&p  +  vs  -  y^gXl^p  +  L^p  +  ^Lg)] 

(30) 


where  the  dashes  indicate  that  for  non-linear  aerodynamics  the  expression  is 
only  valid  for  small 


With  the  gain  parameters  proposed  for  fj  =  J>t  5  and  10  c/s,  given  in 
Tables  V  and  VI  expression  (291  gives 

+  ^  [1  +  0.338e"°,508t  -  3.1l9"3,29t  sin(l.507t  +  25.5°)J  for  f^  =  3% 

*  *  [1  +  0.341  e“°*8l3t  -  3.31e^,25t  sin(2.23t  +  23.9°)]  for  f^  »  5  c/s 

y  a  ^  [1  +  0.324e’"1*37t  -  3.28e~9’25t  sin(3.95t  +  23.8°)]  for  f-j  =  10  o/s. 


These  expressions  are  plotted  in  Pig. 25. 

Por  the  gain  factors  corresponding  to  f  *  a  5  c/s  and  for  the  aerodynamic 
parameters  corresponding  to  zero  and  10°  incidence  at  U  ■  300  ft/seo, 
expression  (30)  gives 
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f  »  [1  +  0.0228e“°‘°57t  +  0.0903e~°*988t  -  4.26e^*53t  + 

at  zero  incidenoe  and  , 

^  [1  +  0.0693Q_0*1^t  -  0.0259s"0, 975t  -  9.36e"5*29t  sin(0.592t  + 6.40] 
at  10°  incidence.  These  two  expressions  are  illustrated  in  Fig.26. 

It  can  be  seen  from  Fig.  25  that  in  the  absence  of  aerodynamic  forces 
the  overshoot  is  approximately  the  same  (about  15$)  for  the  three  servo 
bandwidths  considered.  This  is  in  accordance  with  expectation  since  the  loop 
gains  were  ohosen  to  give  a  60°  phase  margin  in  each  case.  The  higher  speeds 
of  response  obtainable  with  the  larger  loop  gains  is  also  evident. 

Comparison  of  Fig.26  with  the  corresponding  curve  (f*  =  5  o/a)  of  Fig. 25 
indicates  that  the  aerodynamic  foroes  have  the  effect  of  reducing  the  over¬ 
shoot  and  the  speed  of  response.  For  instance,  the  time  at  which  the  heading 
angle  passes  the  demanded  value  is  increased  from  about  0.5  sec  to  about 
0.9  sec  and  the  peak  overshoot  is  diminished  from  15$  to  about  5$*  The  over¬ 
shoot  however  decays  more  Blowly  the  longest  time  constant  being  1.2  sec  in 
the  absence  of  aerodynamics  compared  with  7  seo  for  p  o  10°  and  17*5  sec 
for  zero  incidence. 

5  Conclusions 


This  preliminary  design  study  has  been  made  to  determine  suitable  para¬ 
meter  values  for  two  types  of  heading  control  system  for  RTV1 .  Linear  servo¬ 
mechanism  theory  has  been  used  throughout,  the  analysis  being  based  on  the 
assumption  that  the  short  period  response  of  the  system  and  its  components 
can  be  represented  approximately  by  linear  transfer  functions.  The  maximum 
values  of  the  gain  parameters,  consistent  with  a  well  damped  transient 
response,  are  tabulated  in  section  2  for  the  two  types  of  system  and  for 
three  possible  rooket  motor  servo  bardwidths.  Adequate  damping  is  ensured 
by  choosing  the  values  of  loop  gain  to  give  a  phase  margin  of  60°  in  both  rate 
and  position  feedback  loops. 

Comparison  of  the  responses  of  the  two  types  of  system  to  disturbances 
leads  to  the  conclusion  that  the  system  including  integral  of  error  control 
is  to  be  preferred.  For  example,  it  is  shown  that  with  the  latter  system  a 
rooket  motor  bias  results  in  a  steady  state  component  of  velocity  normal  to 
the  desired  trajeotory,  assumed  straight,  whereas  without  integral  of  error 
a  steady  state  normal  acceleration  component  occurs. 

A  limited  examination  of  the  effects  of  aerodynamic  farces  and  moments 
indicates  that  for  the  conditions  assumed  these  will  have  little  effect  on 
the  response  of  the  proposed  system.  With  an  aerodynamically  unstable  con¬ 
figuration  the  system  beoomes  conditionally  stable,  but  it  is  shown  that  for 
the  integral  of  error  system,  with  the  gain  parameters  proposed  for  a  5  c/s 
servo  bandwidth,  the  gain  margins  in  both  directions  are  large.  With  the 
large  stability  margins  proposed  it  appears  to  be  unnecessary  to  make  provi¬ 
sion  for  varying  the  gain  parameters  in  flight  on  acoount  of  aerodynamic 
forces  or  variations  in  the  vehicle  weight,  moment  of  inertia  eto.  The  com¬ 
bined  effect  of  the  latter  is  to  alter  the  loop  gains  by  about  1  db  between 
tanks  full  and  tanks  empty  conditions. 
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Bar  denotes  Laplace  Transform  (e.g,  5) 

Dot  denotes  differentiation  with  respect  to  time  (e.g. 
Dash  denotes  perturbation  quantity  (e.g.  v' ) 

C  slug  ft^  Missile  lateral  moment  of  inertia 


D1  d b 


D2  db 


e^  volts 

e2  volts 

eT  volts 
3 

fg  cyoles/sec 
f j  cycles/ sec 

2 

g  ft/seo 


Gain  of 


Gain  of 


+  p  —  j  at  frequency  ur 

i  r 

p2g0  L1  +  W3W 


:a2m 


at  frequency  « 


Output  from  heading  gyro  pick-off 

Output  from  rate  gyro  pick-off 

Net  input  to  rooket  motor  servo 

Undamped  natural  frequency  of  rate  gyro 

Undamped  natural  frequency  of  rocket  motor  servo 
(servo  bandwidth) 

Acceleration  due  to  gravity 

frequency  variant  part  of  rate  gyro  transfer  function 

frequency  variant  part  of  rooket  motor  servo  transfer 
function 

frequency  \ariant  part  of  missile  transfer  function  ■ 


volts/degree  Heading  gyro  pick-off  gain  factor 
Kg  volts/degree/sec  Rate  gyro  pick-off  gain  factor  (invariant  with  frequency) 
Kj  degrees/volt  Rocket  motor  servo  gain  factor  (invariant  with  frequency) 

seo~^  frequency  invariant  factor  of  missile  transfer  function 

Kj-  volt3/sec/degree  Integrator  gain  factor 


l  ft 


Distance  between  missile  centre  of  gravity  and  rocket 
motor  pivot  point 

Rats  feedback  loop  gain  (=  K^K^K^) 


seo" 


Position  feedback  loop  gain 


Integral  of  error  loop  gain 


("ft 


-  20  - 


SECRET 


i 

SECRET  | 

t 

Technical  Note  No.  GW  397  j| 

! 

t  ■ 

m  slugs 

mass  of  missile  j| 

1 

N  lb  ft 

Total  moment  in  yaw  plane 

4 

N^  lb  ft/ft/seo 

3N 

Aerodynamio  derivative 

\ 

N  lb  ft/rad/seo 

SN 

Aerodynajnio  derivative  ^ 

Ng  lb  ft/rad 

Motor  control  derivative  ! 

do  i 

rad/seo  /ft/seo 

«,  i 

0  j 

-2 

ng  seo 

H  i 

0  | 

N  [in  eqn.  (28) 
only] 

] 

Number  of  encirclements  of  the  point  (-1  +  jO)  j 

P 

Laplace  Transform  variable 

* 

P 

Number  of  poles  of  feedback  transfer  function  in 
right  hand  half  of  the  p-plane 

r  degrees/sec  or 
rad/seo 

Rate  of  turn  of  missile  (a  *) 

! 

4 

t  seo 

1 

Time  1 

t^  seo 

Time  at  which  rocket  motor  moves  from  limit  j 

T  lb 

Rooket  motor  thrust 

1 

U  ft/seo 

Missile  flight  path  velocity 

V  ft/ seo 

Misaile  side-slip  velocity 

x  ft 

Displacement  in  a  horizontal  direction 

y  ft 

Displacement  in  a  vertical  direction  |' 

*■  ' 

1  lb 

Total  normal  force  in  yaw  plane  i 

I  ,  I  and  - 

v  r  6 

w  ay  ay 

Aerodynamic  and  control  derivatives  —  ,  —  and  -^g 

* 

respectively 

yy  and  y6 

I  Ig 

—  and  —  respectively  j, 

Z 

Number  of  zeros  of  denominator  of  closed  loop  transfer 
function  '  ' 

P  degrees  or 
radians 

Missile  angle  of  incidence  j 

6  degrees  or 
radians 

Rooket  motor  def leotian  with  respeot  to  missile  1 

longitudinal  axis  | 

6^  degrees 

Rocket  motor  deflection  bias  j 
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A  degrees 


Limit  of  rocket  motor  def leotion 


*2 

% 

^  degrees 
^2  degrees 

^  degrees 
^  degrees 

*  degrees 

B 

wq  rad/  sec 
ojc  rad/  sec 
.a>2  rad/seo 


Rate  gyro  damping  faotor 
Servo  damping  faotor 


Phase  of  the  term 


Phase  of  the  term 


at  frequency  wq 


~  W3W2,  ~ 

Li  +  WMV 


at 


frequency  uq 

Heading  angle  of  missile  with  respect  to  vertical  datum 

Demanded  heading  angle  (gyro  datum)  with  respect  to 
vertical  datum 

Heading  angle  error  with  respect  to  gyro  datum  (=  ^  -  ty) 


Break  frequency  of  the  term  ( 1  +  p  —  j  (see  Pig. 9) 

\  K3V 

Angular  frequency  at  which  feedback  gain  is  unity 
(cut-off  frequency) 

2%  f2 


Uj  rad/seo  2%  f^ 
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APPENDIX  I 


Details  of  Trajectory  Calculations 


The  following  calculations  of  approximate  trajectories  are  based  on  the 
assumptions  set  out  in  section  3  and  the  equations  given  in  section  3»  1* 

The  gain  parameters  are  given  in  Tables  IV,.  V  and  VI. 

% 

1.  Deviation  due  to  a  3°  Initial  Heading  Angle  Error 

(a)  Control  system  of  Fig.1  (without  integral  of  error) 

The  control  equation  for  this  case,  from  equation  (l6),  i3 

5  =  •  (1-1 ) 

The  assumed  initial  conditions  are  \Jr  a  3°  and  ijr  a  0  giving  60  « 
3K^Kt  degrees.  It  can  be  seen  from  Table  V  that  for  fj  a  3,  5  anl  10  cycles 
per  second,  3K~|K,  =  3*39°,  8.37°  and  25.92°  respectively.  Thus  the  initial 
value  of  demanded  mo  to  r  deflection,  for  f  x  e  5  and  10  c/s,  exceeds  the 
assumed  5°  limit,  and  so  the  actual  defleotion  will  remain  constant  at  the 
limit  value  for  a  short  period  at  the  beginning  of  .flight.  Referring  to 
equation  (13),  it  is  evident  that  the  missile  angular  motion  during  this 
*  period  is  given  by 


•• 

1^  s 


Tl 

G 


A 


-  V 


(1.2) 


where  A  denotes  the  limit  value  of  motor  deflection.  Henoe 


♦  »  -V*  (*-3) 

t  =  *0  -  i\At2  .  (1.4) 

The  time  at  which  the  demanded  motor  deflection  falls  to  A  may  be 
obtained  by  substituting  expressions  (1.3)  and  (1.4)  in  the  control  equation 
(i.l)  anl  solving  for  t.  Hence 


ID. 

where,  as  before,  hj  =  KgKjK^  and  1^  a  and  only  real  and  positive 
values  of  t^  are  admissible. 

The  trajeotory  during  the  period  0  to  t^  is  obtained  from 


*  «  I  (*  +  A) 


and 


(1.6) 
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whence  the  displacement  (x^ ,y^ )  at  time  t^  is  given  by 


and 


(1.8) 


(1-9) 


The  corresponding  values  of  and  ih  are  obtained  from  expressions 
(l.3)  and  (1.4).  Table  VIII  gives  the  conditions  at  time  for  the  three 
cases  with  the  gain  parameters  specified  in  Table  V. 


Table  VIII 


Conditions  at  time  t^ 


Servo 
bandwidth 
f^  c/s 

*1 

sec 

X1 

ft 

*1 

ft 

*1 

ft/sec 

*1 

,ft/sec 

+1 

degrees 

*1 

deg/sec 

61 

degrees 

3 

0 

0 

0. 

0 

0 

3.0 

0 

3-39 

5 

0.0529 

0.0111 

0.034 5 

0.417 

1.30 

2.85 

-5.87 

5.0 

10 

0.129 

0.0648 

0.205 

0.984 

3.17 

2.08 

-14.30 

5.0 

After  time  h]  the  trajectory  is  obtained  by  solving  equations  (13),  (16), 
(18)  and  (19)  with  the  initial  conditions  given  in  Table  VIII.  Transformation 
of  equation  (13)  gives 


-  \ 


(1.10) 


where,  as  before,  the  bar  notation  indicates  Laplace  Transforms.  Similarly, 
equation  (16)  becomes 


t  =  k1k3V  +  ^(p?-  t|)  . 

Solution  of  equations  (1.10)  and  (i.ll)  gives 

(1.11) 

(p+l.)*,  +  *. 

*  »  T 

p  +  L^p  + 

(1.12) 

,  ,  WlP  +  V*lk”V 

am  0  »  n  * 

\(p  +l1p+l1l£) 

(1.13) 

Prom  equation  (18)  we  have 


p2*  =  f  (^6)  +  x.,  +  px1  (1.14) 

and  substitution  of  expressions  (1.12)  and  (1.13)  in  (1.14)  gives 
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T  Ip(K^+L1L2'{r1+L1*1)  +  K|L1  +  (VLlV^  *1  .  X1 

X  =  21  2  /  2  _  771  1  ~  +  2  +  p 

m  P  K^(p  4-L^pfL^)  P 


Finally,  from  equation  ( i  9 ) ,  we  obtain 


/T  \  l'  yi 

r  "  r8)?  +  r~ 


(1.15) 


(1.16) 


The  missile  motion  after  t^ ,  in  terms  of  ijr,  6,  x  and  y  an  functions 
of  time,  is  obtained  from  the  inverse  transformations  of  expressions  (1.12), 
(1.13),  (l«15)  and  (l.l6).  These  functions  are  given  in  Table  IX  and  Fig. lit- 
shows  the  corresponding  trajectories.  In  this  table  t  denotes  time  after 

V 


Table  IX 


Motion  after  time  t-| 


f3 

ijr  degrees 

£  degrees 

3  c/s 

4.23 e*"3,54t  css (3. 57 K). 782) 

4.81a“3’54t  cos(3.57t*0.776) 

5  c/s 

3.40e"5*65t  css(5.51  t-O.476) 

9.55»“5,65t  cos (5.51 *+1.01 9) 

10  c/s 

2.18e*^*975t  o:s(9.54V0.304) 

18.8a"9'979*  cos  (9.54t+1»305) 

f3 

x  Test 

y  feet 

3  c/s 

0.0102  +  0.834t  - 

12.3*2 

„“3.54t  (0.o-ie2  00s  3.57*  ♦  0.252  sin  3.5?*) 

5  c/s 

O.IO77  +  9.573*  “ 

0.0305  +  1.30*  ♦  12.31? 

e"5.65t  (0.0966  ccs  5.51*  +  0.127  sin  5.51t) 

10  c/s 

0.1558  +  Q.486t  - 

0.205  +  3.17*  ♦  12*3t2 

e"5>.?75*  (0.089  cos  9*54t  ♦  0.0407  sin  9.54t) 

(b)  Oontrol  system  of  Fig.  2  (with  integral  of  error) 

It  can  be  seen  from  Table  VI  that  in  this  oase  also  a  3°  initial  heading 
error  demands  an  initial  motor  deflection  greater  than  5°  f°r  the  gain  para¬ 
meters  corresponding  to  =  5  and  10  cycles  per  second.  Hence  for  these 
two  oases  the  angular  motion  during  the  launching  phase  is  again  given  by 
equation  (1.2),  (l.3)  and  (i.lf) .  In  addition  we  have 

l  t  dt  =  tQt  -  l  K4At3  (1.17) 

assuming  that  the  integrator  output  is  held  at  zero  until  the  instant  of 
launch  (t  =  0).  , 
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Substitution  of  expressions  (1.3),  (1. 4)  and  (X.  1 7)  in  (l 7) >  the  control 
equation  for  this  ca3e,  gives  the  following  equation  for  tj ,  the  time  at 
which  the  demanded  motor  deflection  falls  to  A: 

*2  * 6  v)  ‘ +  6  fe "  w) '  °  '  (I-18> 

K5  K5 

•  L5  =  kJ  L6  s  %  • 

The  required  value  of  t-j  is  given  by  the  smallest  real  and  positive  root 
of  equation  (1.18)  and  the  displacement  at  this  time  is  obtained  from  expres¬ 
sions  (1.8)  and  (1.9).  Table  X  gives  the  conditions  at  t-j  for  the  three 
oases. 


where  L 


'1 


Table  X 


Conditions  at  time  t-j 


Servo 
bandwidth 
f3  c/s 

*1 

sec 

*1 

ft 

y1 

ft 

*1 

ft/sec 

^1 

ft/sec 

*1 

degrees 

*1 

deg/seo 

51 

degrees 

3 

0 

0 

0 

0 

0 

3.0 

0 

2.21 

5 

0.0102 

0.0002 

0.0013 

0.035 

0.251 

2.99 

-1.13 

5.0 

10 

0.1057 

0.019 

0.137 

0.354 

2.60 

2.38 

-11.75 

5-0 

After  time  t.  the  missile  trajectory  is  obtained  from  the  solution  of 
equations  (13),  (l/),  (18)  and  (19)  with  the  initial  conditions  given  in 
Table  X.  Differentiating  and  transforming  equation  (l7)  gives 

P*  -  ft,  «  +  ^(pt-*,)  +  ^(p^+rP^)  (1-19) 


and  solving  (1.1.9)  with  (1.10)  gives 


(p2+4j  P+1J1  Iv, )  +  ^(pfrLj)  -Vi 

¥  ~x  V  2 

P  +  IijP  +  L^P  +  I^Lg 


and 


5  => 


LiLg[p*i +  (^  +  i)]+pVi 


Substitution  of  expressions  (l.20)  aid.  (1.21)  in  (X.1'4)  gives 


X  s= 


♦  t. 


LiL2N 

p(1+i)  +  Li 


(1.20) 


(1.21 ) 


2  3  2 

mp  (pr4L.jP  +L1I^p+L1Lg) 


*1  *1 
T  +  "2 

F  P 


and  y  is  again  given  by  expression  (1.16). 
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The  motion  after  time  t^ ,  in  terms  of  ijr,  5,  x  and  y  as  functions 
of  time,  is  derived  from  the  above  equations.  These  functions  are  given  in 
Table  XI,  in  which,  t  denotes  time  after  t<|,  and  Pig.  15  shows  the  corres¬ 
ponding  trajectories. 


Table  XI 


Motion  after  time  t-j 


h 

C/S 

^  degrees 

6  degrees 

3 

e“3.285t  (4.01  cos  1.507t+8.42  sin  1.5>7t) 

e-3.285t  ^  !  ,507t_  5,02  sin  1.5076) 

-  1 .01  e^- 5081 

+  S,0l2e“°*so8t 

5 

8“5.245t  (k>00  0os  2,2361+8.14  sin  2u236t) 

e-5,2456  (Iw95  ^  2^36t-l3^4  sin  2.23a ) 

+  0#05e"°*81^t 

1c 

Q-9.29t  (2.95  00a  3.956-  3.83  sin  3.95t) 

e"9*29*  (4.90  cos  3.956-27.1  sin  3.956) 

-  0.57a"1 ,57t 

+  o.lOe"1*376 

lls 

x  faet 

y  feet 

3 

3.32  -  3.a2o"°*5)8t 

12.362 

+  a  (tonped  oscillatory  tana 

5 

1.36  -  0.0476  -  UJ13e“0'613t 

0.0013  +  0.251 6  +  12.362 

+  a  cfenped  oscillatory  tats 

10 

0,43  -  0.46t  -  o^Se*1 ,37s 

0,137  +  2.606  +  12.362 

+  a  danped  oscillatory  tern 

The  damped  oscillatory  transient  term  of  x  is  not  given  in  detail  in 
the  above  table  as  it  beooraes  negligible  after  one  second. 

2.  Deviation  due  to  a  1°  Motor  Defleotion  Bias 

(a)  Control  system  of  Pig.  1  (without  integral  of  error) 

In  this  case  it  is  assumed  that  a  constant  bias  exists  in  the  oontrol 
system,  such  that  the  actual  rochet  motor  deflection  is  given  by 

6  =  K^t  +  -  6±  ,  (1.23) 

instead  of  by  equation  (l6).  No  limiting  oocurs  in  this  case  since  no 
initial  heeding  error  is  assumed  and  the  bias  6^  is  assumed  equal  to  1  . 
Transformation  of  equation  (l.23)  gives 

6. 

6  .  K^Kjt  +  KgKjPt  ~  y  U*24) 
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and  from  equation  (l3)  we  have 


2f  V  “ 

P  ♦  =  - 

(1.25) 

whence 

r  'hh  . 

'  p(p2  +  tjP  + 

(1.26) 

and 

.  -p  st 

0  **  2  ' 

p  +  LlP  +  LiL2 

(1.27) 

Transforming  equations 

(l8)  and  (19)  and  substituting  gives 

-  '  T'(\  -  P2)  6i  ' 

X  a  - 5 - £  - 

mp^p  +  L^p  +  LjL^) 

(1.28) 

and 

T  -  mg 

y  =  — * 

mp 

(1.29) 

The  trajectories  obtained  by  transforming  equations  (1.28)  and  (1.29)  into 
functions  of  time  are  shown  in  Pig. 1 6  for .the  three  sets  of  parameters  corres¬ 
ponding  to  fjj  -  3,  5  and  10  cyoles  per  second.  Table  XII  gives  ijq  5,  x 
and  y  as  functions  of  time  for  the  three  cases. 


Table  XII 

Motion  due  to  1°  bias  ' 


f3 

C/S 

ijf  degrees 

5  degrees 

0.88  [1  -  e*"?*5^  (00s  3.57*  +  0.932  sin  3.57*)  ] 

e“3*54*  (0,992  Sin  '3,576  -  cos  3.57*) 

0.36  [1  -  a“5^5t  (00s  5.51*  +  t»025  sin  5.51*)] 

9-5.65*  (1^25  sin  5.51*  -  00s  5.51*) 

10 

0.116  [1  -  (cos  9.54*  +  1*045  sin  9.54*)] 

a-9.?75t  (1.045  sin  9.54*  -  cos  9.54*) 

■_  x  feet  . 

y  feet 

0.435*2  -  0.244*  “  0,0054 

+  a  damped  oscillatory  tan 

'  12.3t2 

0.l765t2  -  0«064t  -  o.olo 

+  a  damped  oselllatory  tern 

• 

Oil 

10 

0.057562  "  0.0l2t  -1.0.0045 

+  a  damped  oscillatory  tara 

12.3*? 
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(b)  Control  system  of  Fig. 2  (with  Integral  of  error) 

It  is  assumed,  in  this  oase  that  the  rooket  motor  deflection  is  given  by 


6  =  K5K3  J  «dt  +  +  K -  6±  (1.30) 

instead  of  by  equation  (17)*  Again  no  limiting  occurs  as  all  initial  con¬ 
ditions  are  assumed  to  be  aero  and  64  is  assumed  equal  to  one  degree. 

The  following  expressions  for  the  Laplace  transforms  of  6  and  x  are 

obtained  from  equations  (13),  (18)  and  (I.30):- 


and 


l 


X 


3  2 

P  +  +  L^LgP  + 

2  * 

-  P  &i 

_  _ 

P  +  +  L^p  +  L_jLg 


"(VP  >  i 

'  '  '  *7  £  '  '  7 

mp  (p  +L^  +L1L2p  +  L1Lg) 


(1.31) 


(1.32) 


(1*33) 


y  is  again  given  by  expression  (I.29). 

The  trajectories  obtained  from  equations  (I.29)  and  (1.33)  ©re  shown 
in  Pig.  17  and  Table  XIII  gives  x  and  y  as  functions  of  time  for  the 
three  oases. 


Table  XIII 


Motion  due  to  1°  bias 


f3 

x  feet 

y  feet 

c/s 

3 

-  8.2  +  3»30t  +  8.4®""°*^°^  +  a  damped  oscillatory  term 

12.3t2 

5 

-  1.29  +  0.83t  +  1.31e~^’^^  +  a  damped  oscillatory  term 

I2.3t2 

10 

-  0.143  +  0. 157t  +  0.137e_1*37t  +  a  damped  oscillatory  term 

12.3t2 

>  7 
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TN,  QW.  3  97. 

FIG.I&2. 


heading 


FIG.I.  BLOCK  DIAGRAM  OF  CONTROL  SYSTEM 
USING  HEADING  ANGLE  ERROR  AND  RATE 

FEEDBACK. 


HEADI NQ  GYROSCOPE 


FIG. 2.  BLOCK  DIAGRAM  OF  CONTROL  SYSTEM 
USING  HEADING  ANGLE  ERROR,  INTEGRAL  OF 
ERROR  AND  RATE  FEEDBACK. 
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TN.  QW-3  9  7. 

FIG  .3. 


FIG. 3.  CONVENTION  FOR  MISSILE  DYNAMICS 
IN  YAW  NEGLECTING  AERODYNAMICS. 

(angles  shown  positive) 
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FIG  .4. 
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FREQUENCY 
CYCLES  PER  SECOND 

FIG.4.  FREQUENCY  RESPONSE  CURVES  FOR  ASSUMED  SERVOS  AND 

RATE  GYRO. 
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FIG  5. 
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FIG. 5.  RESPONSE  OF  FREQUENCY  VARIANT  Ge  G3  G4  FOR  THE  RATE 

FEEDBACK  LOOP 
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FIG.  6 
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RA.D/  SEC. 


CONTROL  SYSTEM  IN  FIG.  2. 
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FIG  12. 


9Q  N 


0-1  1*0  w  RADj  SEC.  10 

FIG.  12.  CLOSED  LOOP  RESPONSE  FOR  THE  OUTER  LOOP  OF  THE 

CONTROL  SYSTEM  IN  FIG.  2. 
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CURVE  I  ©  WITHOUT  INTEGRAL  OF  ERROR 
CURVE  H  0  WITH  INTEGRAL  OF  ERROR 
$5*10  CYCLES  PER  SEC.  IKI  BOTH  CASES 
FREQUENCIES  IN  RADIANS  PER  SECOND 
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FIG  .14, 


f, .  servo  «»ndwidth 


FIG.I4.  MISSILE  TRAJECTORIES  RESULTING 
FROM  A  3°  INITIAL  HHEADING  ANGLE  ERROR 
FOR  THE  CONTROL  SYSTEM  OF  FIG.I. 

(WITHOUT  INTEGRAL  OF  ERROR). 
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FIG. 15. 
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FIG.  15.  MISSILE  TRAJECTORIES  RESULTING 
FROM  A  3°  INITIAL  HEADING  ANGLE  ERROR 
FOR  THE  CONTROL  SYSTEM  OF  FIG. 2. 

(WITH  INTEGRAL  OF  ERROR). 
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FIG. 16 
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FIG.  16.  MISSILE  TRAJECTORIES  RESULTING 
FROM  A  CONSTANT  ROCKET  MOTOR 
DEFLECTION  BIAS  OF  1°,  FOR  THE 
CONTROL  SYSTEM  OF  FIG.I. 

(WITHOUT  INTEGRAL  OF  ERROR). 
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FIG.I7.  MISSILE  TRAJECTORIES  RESULTING 
FROM  A  CONSTANT  ROCKET  MOTOR 
DEFLECTION  BIAS  OF  1°,  FOR  THE 
CONTROL  SYSTEM  OF  FIG.2. 

(WITH  INTEGRAL  OF  ERROR). 
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FIG. 18 
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FIG.  18.  AERODYNAMIC  FREQUENCY  RESPONSE (~K4 G^-  ( AT  $ 

A  =  *'950  0  +  6-84^  AT  ft. 

(•  b-n  o  Kxi  \  ^ '  N  Iv. 


TN.QW.  3  97. 

FIG  10 


FIG.  19.  COMBINED  RESPONSE  OF  SERVO  AND  RATE  GYROSCOPE  (K2  Ga  K3G3) 
(SERVO  BANDWIDTH  OF  5c/*.,  Ki=0  507  AND  K3*l.) 
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FIG.22. 
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FIG.24.  NYQUIST  DIAGRAM  OF  THE  OUTER  LOOP  FOR  THE  CONTROL  SYSTEM 
OF  FIG.2  IN  THE  PRESENCE  OF  AERODYNAMICS. 
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FIG.25.  MISSILE  RESPONSE  TO  A  STEP 
FUNCTION  OF  DEMANDED  HEADING  ANGLE 
tf0  FOR  THE  CONTROL  SYSTEM  OF  FIG.2, 
NEGLECTING  AERODYNAMICS. 
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FIG. 26.  MISSILE  RESPONSE  TO  A  STEP  CHANGE 
OF  DEMANDED  HEADING  ANGLE  ^  FOR  THE 
CONTROL  SYSTEM  OF  FIG.  2.  IN  PRESENCE  OF 

AERODYNAMICS. 
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